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Abstract 

We review the latest calculations of the next-to-leading AS = 1 
effective Hamiltonian, relevant for K — > vrvr transitions. Numerical 
results for the Wilson coefficients are given for different regularization 
schemes. Predictions of e'/e, obtained using different approaches to 
evaluate the relevant hadronic matrix elements, are compared. Given 
the present value of the top mass, rrit = (174 it 17) GeV, all the 
analyses, in spite of the large theoretical uncertainties, indicate that 
the value of e'/e is smaller than 1 x 10~^. 



1 Introduction 

The understanding of mixing and CP-violation in hadronic systems is one of 
the crucial tests of the Standard Model. In the last few years considerable 
theoretical and experimental effort has been invested in this subject. 

On the theoretical side, the complete next-to- leading expressions of the 
relevant effective AS* = 1, AS* = 2, AB = 1 and AB = 2 Hamiltonians have 
been computed thus reducing the theoretical uncertainties^ More- 

over, there is now increasing theoretical evidence that the value of the pseu- 
doscalar i?-meson decay constant is large, fs ~ 200 MeV, and that the B B 



^Indeed only the top contribution to the AS — 2 Hamiltonian is fully known, at the 
next-to-leading order. To our knowledge, the other contributions have been only partially 
computed M. 
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parameter Bb is quite close to one. This strongly constrains the Cabibbo- 
Kobayashi-Maskawa parameters and it has remarkable consequences on CP- 
violation in B decays, see refs. ^. Still, the evaluation of hadronic matrix 
elements is subject to large uncertainties, that are particularly severe for 
e'/e, where important cancellations of different contributions occur for large 
values of the top mass. Indeed, a significative reduction of the theoretical 
uncertainty on e'/e would require a substantial improvement in the calcula- 
tion of the hadronic matrix elements, either from lattice simulations or from 
other non-perturbative techniques. 

On the experimental side, more accurate measurements of the mixing 
angles are now available and the mass of the top quark, experimental evidence 
of which has recently been found by CDF |rD||, is constrained within tight 



limits [11|. Still, in spite of very accurate measurements, the experimental 



results for the CP violating parameter e'/e are far from conclusive [|12|, |T3l. A 
better accuracy, at the level of 1x10""^, should be achieved by the experiments 
of the next generation. 

In the following, we briefly introduce the AS = 1 effective Hamiltonian 
and summarize the main results of the next-to-leading calculation of the 
relevant Wilson coefficients. An updated analysis of e'/e, along the lines 
followed in refs. 0-0, is presented. Particular emphasis is devoted to a 
realistic evaluation of the uncertainties. We also compare the results of refs. 
[0]^0 next-to-leading order analysis of ref. |]14|. In section ^, the basic 



formulae, which define the CP-violation parameters e and e', are presented; 
in section 0, the definition of the Cabibbo-Kobayashi-Maskawa matrix and 
the notation used in this work are introduced; in sections | and ^, we give 
several details about the AS = 1 effective Hamiltonian relevant for direct 
CP-violation. In particular, the Wilson coefficients in different regulariza- 
tion schemes are reported. In section 0, we give the formulae which has been 
used to obtain the theoretical predictions; in section ||, the theoretical pre- 
dictions from the more recent analyses are given. Further details, including 
a theoretical discussion of the matching conditions, of the 5-parameters and 
of the uncertainties coming from the choice of Aqqd, the renormalization 
scale, etc. can be found in ref. M. 
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2 CP-violation in K ^ nn decays 

In this section, we introduce the parameters e and e' that describe CP- 
violation in the neutral kaon-system. In the following, we assume CPT 
symmetry. A comprehensive discussion of the general case, including CPT- 



violation, can be found in ref. [|T5|. 



There are two possible sources of CP- violation in the decays of the neutral 
kaons into two pions. CP-violation can take place both in the kaon mixing 
matrix and at the decay vertices. Let us consider the mixing first. The most 
general CPT-conserving Hamiltonian of the K^-K^ system at rest can be 
written as 

where the bra (ket) can be represented as the two component vector {K^\ = 
(1,0), {K^\ = (0,1), M is the "mass" matrix and F is the "width" matrix. 
Both Mo and Fq are real. 

Notice that there is some freedom in the definition of the phases of the 
kaon field. In particular, one can make the change 

e*"|ir°) , ^ e'^'lK^) . (2) 

Correspondingly, the off-diagonal matrix elements of any operator X, acting 
on the K^-K^ system, undergo the changes 

Xi2 ^ e-2-Xi2 , X21 ^ e2-X2i . (3) 

This arbitrariness enters in some popular definitions of the CP-violation pa- 
rameters. For definiteness, we choose a particular phase convention, namely 
we require that the CP operator is given by 



In this case, 

\K±) = ^ {\K') ± \K')) (5) 
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are the CP eigenstates. In the presence of CP- violation, [CP, i/] 7^ and the 
eigenstates of the Hamiltonian H are not CP eigenstates. We introduce the 
parameter e which defines the eigenstates of H as 

' + |ep \/l + |eP 



The corresponding (complex) eigenvalues are denoted as 

1 



\l = mL-iTL/2 = -{Hi^ + H22)-\jHi2H2i 



(7) 



\s = ms-tTs/2 = -{Hn + H22) + ^/Hi2H2i. 



In the phase convention @, the parameter e controls the amount of CP- 
violation, namely the CP symmetric limit is recovered for e ^ 0. We can 
explicitly write e in terms of the matrix elements of H 

H21 — H12 
AA — H12 — H21 

where AA = A^ — A5. 

Experimentally CP-violation is a small effect, i.e. |e| <^ 1. For this 
reason, one can simplify eq. (|]) to obtain 

iImMi2 + Imri2/2 

f ~ 

AA 

with 



(9) 



AA = AM-iAr/2 
AM = -2ReMi2 (10) 
AF = -2ReFi2. 



Moreover, since AM/AF = -0.9565 ± 0.0051 ^ -1, eq. (|) becomes 

_ ^l+i ImMi2 1 - i ImFi2 
^ ~ 2 2ReMi2 2~2Ref^ ' 



(11) 
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In view of the following discussion of the CP-violation parameters, let us 
introduce amplitudes of the weak decays of kaons into two pions states with 
definite isospin 

Aje''' = {nn{I)\Hw\K^) , (12) 

where / = 0, 2 is the isospin of the final two-pion state and the 5/'s are the 
strong phases induced by final-state interaction. Watson's theorem ensures 
that 

A*je''' = {7nT{I)\Hw\K') . (13) 

Direct CP-violation, occurring at the decay vertices, appears as a difference 
between the amplitudes {ttttIH^^IK^) and {7m\H\Y\K^) ■ This corresponds to 
a phase difference between Aq and 

One introduces the parameter e' to account for direct CP-violation. A 
convenient definition is 

e' - 



{mr{0)\Hw 


Ks){nni2)\Hw 


Kl) - {mr{0)\Hw 


KL){7Tn{2)\Hw 


Ks) 




V2{ 


n'K{0)\Hw 


Ks? 



where u = Reyl2/Reylo. Equation (|1J) is obtained in the approximation 



IihAq -C ReAo, Im742 -C Re742 and also u; <^ 1, as a consequence of the 
AI = 1/2 enhancement in kaon decays; e' is independent of the kaon phase 
convention. On the contrary, the parameter e, defined in eq. (^), depends 
on the choice of the phase. Under a redefinition of the phases as in eq. (|^), 
e changes as 

—2 sin a -|- ecosa , 

:—. 15 

cos a — le sm a 

and the CP-symmetric limit does not correspond to e ^ 00. 

Another parameter, which is independent of the phase convention and 
accounts for CP-violation in the mixing matrix, can be defined in terms of 
the K TTTT transition amplitudes 

{mT{0)\Hw\KL) i sin 00 + e cos 00 .ImAo 

e = — j = ~ e -|- I — , (ibj 

{■ktt{0)\Hw\Ks) cos 00 -|- ie sin 00 Re^o 



^In this case, the states \K±) are not CP eigenstates. 
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where Aq = \Ao\e'^^° and the last expression is obtained in the approximation 
e, 00 ^ 1- The two definitions, eqs. (H) and (|16D, coincide in the Wu-Yang 
phase convention, ImAo = 0. One can check that 0o changes with the phase 
convention as 0o ^ 0o + and that e is invariant. 

From unitarity, one has 

T,2 = J2'^''^iMK - En){K'^\Hw\n){n\Hw\K^) . (17) 

n 

Given the dominance of — > vrvr (0) decay, one obtains the relation = 
(Aq)"^. From eqs. (|lT]) and ([T6|), one has 

where ^ = — ImAo/ReAo. In the Cabibbo-Kobayashi-Maskawa phase con- 
vention, the ^ contribution is small and can be safely neglected. 

To make contact with the experiments, one defines the two amplitude 
ratios 





Hw 


Kl) 




Hw 


Ks) 



(tt+tt 


Hw 


Kl) 


(tt+tt" 


Hw 


Ks) 



Neglecting small terms, one has 

r^oo - e - 2e' , 77+. ~ e + e' , (20) 

namely 



-(7) -K'-S)- '^^^ 

Expressing r/oo and 7]^ in terms of the corresponding widths 





ttOttO) 


T{Ks- 


^ ttOttO) 




TT^TT^) 


T{Ks- 


TT+TT^) 



^+-r = (22) 



eq. (|2l|) gives the CP-violation parameters in terms of measurable quantities. 
Notice that e'/e is approximately real, since experimentally 62 — So ^ —tt/4. 
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3 The CKM matrix 



In the Standard Model, the basic quark charged- current interactions are de- 
scribed by the Lagrangian 
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quark— W 



uaf.{i - i^WijdjW + h.c. 



(23) 



where Ui are the charged 2/3 quarks {u, c, t), dj the charged —1/3 quarks 
{d, s, h) and g is the SU{2)i weak couphng constant [Gpj^ = g^/SM\ 



where Gp is the Fermi constant). V is the unitary CKM matrix |jT6 
useful parametrization is O, |18i 



A 



V 



us 

Vcs 



'ts Vtb 
CeCcr 
— SgCr — CgS^S-re^^ 



(24) 



SgCa 

— CgSr — C-rSgSo-e^^ 




In eq. (p^ , 6, a and r are quark mixing angles (in particular, 6 corresponds 
approximately to the Cabibbo angle); Cg, Sg, etc., mean cos6', sin^, etc.; S 
is the CP- violating phase. Experimental determinations of \Vud\, \Vcb\ and 
\Vub\ from K and B decays show that there is a hierarchy in the mixing 
angles, so that the CKM matrix can be empirically expanded in powers of 
A = 5*61 ^ 0.22 ||19|. Up to and including terms of order A^ (A^) for the real 
(imaginary) part, V is given by 



V 



1-^ 

2 

2x4/ 



A 



1 



2 



AX' 



1- [1-t)(p + 'V) -AX^l + X\p + IT])) 



AX' (p - if]) \ 
AX^ 
1 

(25) 

where Sr = AX^ and S'^-e = AX'{p — irf). In this particular (quark) phase 
convention, the imaginary part of the matrix appears at order A^. 

The unitarity of the CKM matrix implies 



(26) 
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Figure 1: The unitarity triangle in the plane. 
In particular, considering the condition 

KT^Kd + KtKd + v*v,, = 

in the approximation Vud — — 1, one obtains 



+ 



1 = 0. 



(27) 



(28) 



This relation identifies a triangle in the p-r] plane (see fig. p. The angles of 
this triangle, a, (3 and 5, are measures of CP- violation. 

Recent phenomenological analyses of the CKM matrix elements can be 
found in refs. P, Eol A brief discussion of these analyses together with 



the numerical results, can be found in section R. 



4 The bare AS* = 1 effective Hamiltonian 



Weak decays of light hadrons are more conveniently studied using the Wil- 
son operator product expansion (OPE) [E^. With the OPE, it is possible to 
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introduce an effective Hamiltonian, written in terms of renormalized local op- 



erators and of the corresponding Wilson coefficients |^-||26[. Short-distance 
strong-interaction effects are contained in the coefficients and can be com- 
puted in perturbation theory, because of asymptotic freedom. Long-distance 
strong-interaction effects are included in the hadronic matrix elements of the 
local operators and must be evaluated with some non-perturbative technique 
(lattice, QCD sum rules, etc.). The convenience of the effective Hamiltonian 
approach is that all known non-perturbative methods are usually able to 
predict matrix elements of local operators only. In this section we introduce 
the bare AS = 1 effective Hamiltonian, the renormalization of which will be 
discussed in the next section. 

At second order in the weak coupling constant and at zero order in the 
strong coupling constant, the AS = 1 effective Hamiltonian can be written 
in terms of a local product of two charged currents 



Gp 

T-Cw^^ = 7k l-^u i^L'Jf.UL) (uLj^dL) + Ac 



V2 



u 



72 



Au {SaUa)(^V-A) (W/?^^/?) (V-A) + >^c {u ^ c)\ , (29) 



where (sa(io)(y_^) = ■Sq7^(1 — 75)^0; a and /? are colour indices and the sum 
over repeated indices is understood. We have introduced the notation 

K = y,<iv;, (30) 

for q = u,c, t. In terms of the Ag, the unitarity condition of the CKM matrix 
can be written as 

A„ + Ac + At = . (31) 

Equation (|29|) has been obtained from the original theory, by neglecting 
all masses and momenta with respect to Mw- In practice, the effective 
Hamiltonian is obtained by taking 1/(M^ — q^) — >• 1/M^ in the T-product 
of the two charged currents and by putting the m, d and s masses to zero. In 
order to discuss CP-violation, it is convenient to write Ti^f^^ as 

n%'=' = -Xu^ [(1 - r) (g^ - g^) + rg^] , (32) 

where r = —\t/\u contains the CP-violating phase and 
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5 QCD corrections 



Strong interactions play a crucial role in non-leptonic weak decays. The per- 
turbative short-distance effects, included in the calculation of the Wilson co- 
efficients, may be very important because of the presence of large logarithms 
~ Q;"ln™(Mvi///i), where /i is a scale of the order of the mass of the decaying 
hadron. For an accurate estimate of the short-distance contributions, the 
large logarithms have to be resummed to all orders using renormalization 
group (RG) techniques. 

The starting point is the T-product of the two weak currents expanded 
at short distances in terms of local operators. Taking into account the renor- 
malization effects due to strong interactions, we write 

2 

{F\n^'=^\I) = ^-j d'xD^^ {x,Mw) (F|T(j^(a;), J,(0))|/) = 

-^Ea(/i)(F|g,(/.)|/) + ..., (34) 

where {F\ and |/) are the generic final and initial states; the Qiin) form a 
complete basis of operators renormalized at the scale /i; the Cj(/i) are the 
corresponding Wilson coefficients and the dots represent terms which are sup- 
pressed with respect to the dominant ones as powers of Aqqj) /M^r {ml/M'^r 
for S-decays). The effective Hamiltonian is independent of renormalization 
scale /i. On the lattice, the renormalization scale can be replaced by the 
inverse lattice spacing and the effective Hamiltonian can be expressed in 
terms of bare lattice operators . The OPE in eq. (0) must be valid for all 
possible initial and final states. This implies that the effective Hamiltonian 
is defined from an operator relation 

n%'=' = E a(^)g.(^) = -^Q^i^^) ■ cif^) ■ (35) 

The important features of 'H^P^^ are the following: 

• the Wilson coefficients can be calculated using (RG-improved) per- 
turbation theory, provided that one chooses a sufficiently large renor- 
malization scale fi ~ 2-3 GeV ^ Aqcd- In the leading logarithmic 
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approximation (LLA), all terms of O («s(/i)"log(Mvi//yu)") are taken 
into account; 

• all non-perturbative effects are contained in the matrix elements of the 
local operators, the calculation of which requires a non-perturbative 
technique. 

Since TC^^^, eq. (pSD , is independent of fi, the coefficients C(/i) = 
(Ci(/i), C2(/i), . . .) must satisfy the RG equations 

^'"^2Ci^^) = l'f^ci^^) , (36) 



which can be more conveniently written as 



^',^ + /3K)-^- Vk)) (?(,<) = 0, (37) 



where 

f3{as) = (38) 

is the QCD /5-function and 

7(a.) = (39) 

is the anomalous-dimension matrix of the renormalized operators; Z is de- 
fined by the relation which connects the bare operators to the renormalized 
ones, Q{j2) = Z~^{jj, as)Q^ . 

The solution of the system of linear equations ( P7D is found by introducing 
a suitable evolution matrix ?7(/x, M^) and by imposing an appropriate set of 
initial conditions, usually called matching conditions. The coefficients C(/i) 
are given by^ 

Cii^) = Uilj,,Mw)CiMw), (40) 

with 

U {nil, 1712) =Ta,exp[ ^5r^7^("^)h (^1) 



^ The problem of the thresholds due to the presence of heavy quarks with a mass 
Mw 3> niQ Aqcd will be discussed below. 



11 



is the ordered product with increasing couphngs from right to left. The 
matching conditions are found by imposing that, at /i = Mw-, the matrix 
elements of the original T-product of the currents coincide, up to terms 
suppressed as inverse powers of Mp^/, with the corresponding matrix elements 
of Ti^^^. To this end, we introduce the vector T defined by the relation 

(a|T(jtj)|/5) = _^(«|Q^|/3)o . f + . . . (42) 

where {a\Q'^\(3)Q are the matrix elements of the operators at tree level. We 
also introduce the matrix M(/i) such that 

{a\n%'='\(3) = -%(«|g^(;x)|/5)C(/x) = 



-^(«|g^|/5)oM^(/i)C7(/i). (43) 

In terms of T and M, the matching condition 

{a\T{J^jm = {a\n^'='\P) (44) 

fixes the value of the Wilson coefficients at the scale M\y 

C{Mw) = [M^{Mw)]~^f . (45) 

Notice that the matching could be imposed at any scale /i, such that large 
logarithms do not appear in the calculation of the Wilson coefficients at the 
scale /i, i.e. aslnMvi///i ^ 1- 

Equation (^Oj) is correct if no threshold corresponding to a quark mass 
between and Mw is present. Indeed, as a^, 7 and [3{as) all depend on 
the number of active flavours, it is necessary to change the evolution matrix 
U deflned in eq. (^T]), when passing the threshold. The general case then 



corresponds to a sequence of effective theories with a decreasing number 
of "active" flavours. By "active" flavour, we mean a dynamical massless 
(yU ^ mg) quark fleld. The theory with k "active" flavours is matched to the 
one with k + 1 "active" flavours at the threshold. This procedure changes 
the solution for the Wilson coefficients. For instance, if one starts with five 
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"active" flavours at the scale My/ and chooses nic <^ <^ m^, the Wilson 
coefficients become 

= W[^i, Mw]C{Mw) = Mfi, mh)U5imb, Mw)C{Mw) • (46) 
The inclusion of the charm threshold proceeds along the same lines. 

6 The operators of the AS* = 1 effective Hamil- 
tonian 

So far, we have presented the exact solutions of the renormalization group 
equations for the Wilson coefficients. In practice, it is only possible to calcu- 
late the relevant functions in perturbation theory. For illustrative purposes, 
we consider the calculation of the AS = 1 effective Hamiltonian at the leading 
order in QCD. The bare Hamiltonian is given in eq. (|29|). In the presence 
of QCD interactions, other operators appear in the Wilson expansion. A 
complete basis is given by the following operators 



Qi = 




(V- 


-A)iui3Ui3)(V-A) 


Q2 = 


{Sad/3, 




-A)iUl3Ua){V-A) 


Qs = 


i^Sada, 


iv~ 


-A) '^i^l3(l/3)iV-A) 

q 


Q4 = 


{Sad/3, 


(V- 


-A) '^{Q/3(la)iV-A) 

q 


Q5 = 




iv- 


-A) '^{Qp<lp)iV+A) 
q 


Qe = 


[Sadjs] 


(V- 


-A)^iqf3qa)(V+A) 

q 


Ql = 


(ySada^ 


iv- 


-A)iCf3Cp)(V-A) 


Ql = 


{Sadp, 


iv- 


-A)iCl3Ca){V-A) ■ 



The q index runs over the "active" flavours. The above operators are gener- 
ated by gluon exchanges in the Feynman diagrams of fig. |^. In particular, Qi 
is generated by current-current diagrams and Qs-Qe are generated by pen- 
guin diagrams. The choice of the operator basis in not unique, and different 
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possibilities have been considered in the hterature If the electromag- 

netic correction, are also taken into account, the operator basis enlarges to 
include the following operators 





3 

2' 


[Sada){V-A) (^qiWp) (V+A) 

q 


Qs 


3 
2' 


[Sadi3){V~A) (^qi.Wa){V+A) 

q 


Q9 


3 

2' 


[Sada){V-A) ^ (iq{.Wp) {V -A) 

q 


Qio 


3 

2' 


[Sadfi){V-A) X! '^q{Wo){V-A) 

q 


Below the bottom threshold. 


the following relation holds 




Qio 


- - <54 + = , 



so that there are nine independent operators. The basis is further reduced 
below the charm threshold by using the relations 

Q4-Q3-Q2 + Q1 = 

Qg-^gi + igg = 0. (50) 



All the operators considered above are dimension-six operators. In prin- 
ciple, two dimension-five operators 

Vll — ^Q^2^^^'^^{V—A)'^ t^^ 

Q12 = Y^msSa(T'^y_A)ti(^dpG^^ (51) 

should also be included in the operator basis. The matrix elements of Qn 
and Q12, however, enter only at 0{p'^) in chiral perturbation theory. Since the 
phenomenological analysis presented in the following is only valid up to terms 
of 0{p^), we do not need to include the contribution of the dimension-five 
operators in the calculation of e'/e. The effect of these operators on e'/e has 
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Figure 2: One-loop corrections to the AS=1 effective Hamiltonian. 



recently been analysed in ref. |2^. Other operators of lower dimensionality 
(e.g. two-fermion operators) are also potentially present. However, it can 
be shown that their effect can be reabsorbed in a suitable redefinition of the 
fermion fields and by diagonalizing the quark mass matrix at first order in 

Gf H-H- 



In summary, the AS = 1 effective Hamiltonian, renormalized at a scale 
rric, can be written as 



+rQ(/x)^C7(/i)} 

where, in order to find the Wilson coefficients to a given order in a^, "we have 
to calculate eqs. (^), (^) in perturbation theory. 

The exphcit expressions of 7s (a^) and in the LLA 

2 



(52) 



An 



7: 



(0) 



Air 



(3o 



(53) 
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can be found for example in ref. p. In eq. (PD, using 7^°^ and Po, one 
obtains 

^"'■^'-)-(^) • 

At this order, the matching conditions are triviah M, eq. (^), is the 
identity matrix; T, eq. (^2|), has all vanishing components with the only 
exception of T2 = 1. Thus the Wilson coefficients at the leading order for 

rric ^ ^ rrib are given by 

y as' (/i) / \ as' [mb) J 

(55) 

with C2{Mw) = 1 and all the other Wilson coefficients at the scale Mw 
vanish. 

In the next-to-leading logarithmic approximation (NLLA), one proceeds 
along the general scheme described above. In this case, all quantities entering 
in the matching procedure have to be computed at order as (oe for the 
electromagnetic case). The /3-function and the anomalous dimension matrix 
have to be computed at second order in the coupling constants. Thus, for 
example, the anomalous dimension matrix in the NLLA has the form 

where 0{al) corrections have been neglected. We will not give here any de- 
tails of the NLLA calculations. They can be found in refs. At the 
next-to-leading order, it is necessary to solve numerically eq. (pTf). Table 
1^ contains the coefficients, calculated at the leading (LO) and at the next- 
to-leading (NLO) order, using the 't Hooft-Veltman (HV) and the naive 
dimensional (NDR) regularization schemes, for different values of the renor- 
malization scale fi. The errors in the table take into account the varia- 
tion of the values of the coefficients due to Aq^^, = (330 ± 100) MeV and 
nit = (174 ± 17) GeV. Notice that the next-to-leading Wilson coefficients 
and operators both depend on the regularization scheme, while the effective 
Hamiltonian is scheme- independent up to terms 0{al). Actually the depen- 
dence of the effective Hamiltonian on the regularization scheme, due to the 
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unknown next-to-next-to-leading terms, can be estimated and contributes to 
the uncertainties in the prediction of e'/e, see ref. §]. 

The coefficients in table |I] have been computed independently by the 
Munich group |14|. The definition of the renormalized operators in the 

This is due to 
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HV scheme used here differ from those defined in ref. 
the different way of taking into account the two-loop anomalous dimension 
of the weak current, which does not vanish in the HV calculation. One can 
relate the HV coefficients of table |1] (C) and those of ref. (C)- The 
relation is 



(57) 



where 



7j 



2Nr 



-/3o 



(58) 



Once these differences in the definition of the renormalized operators and the 
reduction of the operator basis, eq. (|49|), are properly taken into account, 
the numerical results presented here agree with those of ref. |T4 



7 Relevant formulae 

In order to estimate e'/e, we have to constrain the CP- violating phase 5 in 
the CKM matrix, by using the available experimental information. To this 
end, we consider the CP-violating term in the K^-K^ mixing amplitude and 
the CP-conserving term for B^-B^ mixing. In the following, we present all 
the formulae used in our analysis, namely the expressions of e, Xd and e'/e 
from the AS' = 2, AB = 2, AS'=1 effective Hamiltonian, respectively. 

The effective Hamiltonian governing the AS* = 2 amplitude is given by 

= ^.M^isj^^df {xlF{x,) + X]F{xt) + 2A,A*F(x„x,)} , (59) 

where Xq = m'^/M^r and the functions F{xi) and F{xi,Xj) are the so-called 
Inami-Lim functions [^, including QCD corrections 0; F{xt) is known at 
the next-to-leading order and has been included in our calculation. From 
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LO 


NLO HV 


NLO NDR 




At = 1.5 GeV 


Ca 
Cs 
Ci 
Cs 
C(i 
Ct 
Cs 
Cg 


(-4.22 ± 0.65 ± 0.00) x 10"^ 
(11.62 ±0.38 ±0.00) X 10-1 

(1.99 ±0.35 ±0.00) X 10-2 
(-4.16 ±0.56 ±0.02) X 10-^ 

(1.19 ±0.12 ±0.00) X 10-2 
(-0.66 ±0.13 ±0.00) X 10-1 

(0.16 ±0.04 ±0.19) X 10-3 

(0.63 ±0.14 ±0.16) X 10-3 
(-6.77 ±0.27 ±0.71) x 10"^ 


(-3.91 ±0.51 ±0.00) X 10-^ 
(106.13 ± 0.82 ± 0.00) X IO-2 

(2.17 ±0.41 ±0.00) X 10-2 
(-4.51 ±0.60 ±0.01) X 10-2 

(1.37 ±0.15 ±0.00) X 10-2 
(-0.63 ±0.11 ±0.00) X 10-1 
(-0.04 ± 0.00 ± 0.17) X 10-3 

(0.97 ±0.19 ±0.15) X 10-3 
(-6.32 ± 0.37 ± 0.64) x 10-3 


(-3.80 ±0.55 ±0.00) X IQ-i 
(11.95 ±0.35 ±0.00) X 10-1 
(2.60 ±0.52 ±0.00) X IO-2 
(-0.63 ±0.11 ±0.00) X 10-1 
(10.52 ±0.61 ±0.01) X 10-3 
(-0.93 ±0.21 ±0.00) X 10-1 
(0.02 ±0.06 ±0.20) X 10-3 
(1.06 ±0.26 ±0.19) X 10-3 
(-7.24 ±0.19 ±0.73) x 10-3 




^l = 2 GeV 


Ci 
Ca 
C3 
C4 
Cs 
Ce 
C7 
Cg 

Cg 


(-3.47 ±0.44 ±0.00) x 10"^ 
(11.16 ±0.23 ±0.00) X 10-1 
(1.59 ±0.23 ±0.00) X 10-2 
(-3.50 ±0.40 ±0.01) X 10-2 
(10.40 ± 0.94 ± 0.04) X 10-^ 
(-5.23 ± 0.80 ± 0.03) x IO-2 
(0.18 ±0.02 ±0.19) X 10-3 
(0.50 ±0.08 ±0.12) X 10-3 
(-7.03 ± 0.20 ± 0.74) x 10-3 


(-3.29 ±0.37 ±0.00) x lO-i 
(104.13 ± 0.54 ± 0.00) X IO-2 

(1.73 ±0.26 ±0.00) X 10-2 
(-3.82 ± 0.44 ± 0.01) X IO-2 

(1.20 ±0.11 ±0.00) X 10-2 
(-5.08 ± 0.72 ± 0.03) x IO-2 

(0.01 ±0.00 ±0.18) x 10-3 

(0.77 ±0.12 ±0.12) x 10-3 
(-6.71 ± 0.27 ± 0.68) x 10-3 


(-3.13 ±0.39 ±0.00) X 10-1 
(11.54 ±0.23 ±0.00) X 10-1 
(2.07 ±0.33 ±0.00) X IO-2 
(-5.19 ±0.71 ±0.01) X 10-2 
(10.54 ±0.16 ±0.02) X 10-3 
(-0.72 ±0.13 ±0.00) X 10-1 
(0.01 ±0.04 ±0.20) X 10-3 
(0.81 ±0.16 ± 0.14) X 10-3 
(-7.49 ±0.15 ±0.75) x 10-3 




^1 = 3 GeV 


Ci 
Ca 
C3 
C4 
Cs 
Ce 
C7 
Cg 
Cg 


(-2.68 ± 0.28 ± 0.00) x 10-^ 
(10.71 ±0.12 ±0.00) X 10-1 

(1.20 ±0.14 ±0.01) X 10-2 
(-2.78 ±0.26 ±0.01) X IO-2 

(8.60 ±0.67 ±0.04) x 10-3 
(-3.89 ± 0.47 ± 0.03) x IO-2 

(0.22 ±0.01 ±0.20) X 10-3 

(3.73 ±0.50 ±0.92) x 10"'' 
(-7.29 ±0.14 ±0.77) x 10-3 


(-2.59 ± 0.25 ± 0.00) x lO-i 
(101.88 ± 0.24 ± 0.00) X IO-2 

(1.29 ±0.16 ±0.00) X 10-2 
(-3.06 ± 0.29 ± 0.01) X 10-2 

(9.99 ±0.84 ±0.04) x 10-3 
(-3.84 ± 0.44 ± 0.03) x IO-2 

(0.08 ±0.00 ±0.19) X 10-3 

(5.84 ±0.74 ±0.91) x lO"* 
(-7.04 ± 0.18 ± 0.73) X 10-3 


(-2.41 ±0.25 ±0.00) X 10-1 
(11.12 ±0.14 ±0.00) X 10-1 

(1.56 ±0.19 ±0.01) X 10-2 
(-4.10 ±0.46 ±0.00) X 10-2 

(9.75 ±0.46 ±0.03) x 10-3 
(-5.33 ± 0.73 ± 0.03) x IO-2 

(0.04 ±0.02 ±0.20) X 10-3 

(0.59 ±0.09 ±0.10) X 10-3 
(-7.66 ±0.11 ±0.78) X 10-3 




fi = 4.5 GeV 


Ci 
Ca 
C3 
C4 
Cs 
Ce 
C7 
Cg 
Cg 
Cio 


(-2.07 ±0.19 ±0.00) X 10-^ 
(103.85 ± 0.67 ± 0.00) x IO-2 

(1.13 ±0.11 ±0.03) X 10-2 
(-2.42 ± 0.20 ± 0.02) X 10-2 

(7.06 ± 0.49 ± 0.04) x 10-3 
(-2.94 ± 0.29 ± 0.02) X IO-2 

(0.29 ±0.01 ±0.20) X 10-3 

(2.88 ±0.32 ±0.69) x 10-"' 
(-0.96 ±0.01 ±0.10) X 10-2 

(2.12 ±0.19 ±0.24) X 10-3 


(-2.03 ± 0.18 ± 0.00) X 10-1 
(1001.15 ± 0.43 ± 0.00) X 10-3 

(1.19 ±0.12 ±0.03) X 10-2 
(-2.67 ±0.23 ±0.02) x IO-2 

(8.36 ±0.66 ±0.04) x 10-3 
(-2.97 ± 0.29 ± 0.02) x IO-2 

(0.17 ±0.00 ±0.19) x 10-3 

(4.53 ± 0.49 ± 0.70) x lO-* 
(-0.93 ± 0.00 ± 0.10) x 10-2 

(1.95 ±0.17 ±0.24) x 10-3 


(-1.84 ±0.17 ±0.00) X 10-1 
(108.20 ± 0.87 ± 0.00) x IO-2 

(1.40 ±0.14 ±0.03) X 10-2 
(-3.51 ±0.34 ±0.02) X IO-2 

(8.75 ±0.52 ±0.04) x 10-3 
(-4.09 ± 0.46 ± 0.02) x IO-2 

(0.09 ±0.01 ±0.20) X 10-3 

(4.39 ±0.52 ±0.77) x 10-"' 
(-0.97 ±0.01 ±0.10) X 10-2 

(1.86 ±0.16 ±0.23) X 10-3 



Table 1: Wilson coefficients of the AS = 1 effective Hamiltonian at ji = 
1.5,2,3,4.5 GeV. For fi < mi,, the relation has been used to reduce the 
operator basis. We take A'^qIij = (330 ±100) MeVandnit = (174 ±17) GeV. 
The values of the coefficients shown here correspond to the central values of 
these parameters. The first error is due to the uncertainty on Aqcd, the 
second is due to nit. 
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eqs. ([T8| ) and (^91), one can derive the CP-violation parameter 

\e\^=o = C^BkA"^ sin S {F{xc,xt) + 

F{xt)[A^\\l-acos6)]-F{x,)} , (60) 

where 

GlfiM^ (61) 

In eq. (|60D, p = crcos5, = (Tsin(5 and A, A, p and 77 are the parameters of 
the CKM matrix in the Wolfenstein parametrization [T^. Bk is the renor- 
mahzation group invariant S-factor, the definition of which at the leading 
order is 

{K\ (s^^df \K) = yiMla,{pfBK • (62) 
The AS = 2 effective Hamiltonian is given by 



2 



'^efT = ^,Ml,\l{ln1dfF{x,) . (63) 
Here At = VtdV^l. From eq. (^), one finds the B^-B^ mixing parameter 
X, = ^ = Cs^BsA'xHl + a' -2acos6)F{x,), 

i IVI B 

- 6vr2 ' ^^^^ 

where Bb is the 5-parameter relevant ioi B — B mixing, the definition of 
which is analogous to the Bk one. 

We can write e' as 



V2 ReAo 



'uj'\lmA2y - (1 - Qib) ImAol . (65) 



With respect to eq. ([l4|), we have here explicitly written the isospin-breaking 
contribution fliB, see for example ref. pO|, 



(ImAa)' = (ImAa) - ^jBiujlmAo) . (66) 
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To compute ImAo and (1111^2)', we need the hadronic matrix elements of 
the operators Qi between a kaon and two pions. Usually they are given in 
terms of the so-called S-parameters: 

{TTTTil = 0)\Q,{fi)\K) = By\i2){7T7r{I = 0)\Q,\K)viA 

(7r7r(/ = 2)|g,(/.)|/r) = B'f{i,){7r7r{I = 2)\Qi\K)viA, (67) 

where the subscripts VIA means that the matrix elements are computed in 
the vacuum insertion approximation. The relevant VIA matrix elements can 
be expressed in terms of three quantities 

X = U{MI-Ml), (68) 

\ms{iJ,) + nidiiJ,) J V ''T^sW J 

Z = 4||y^-ljF. (70) 

Prom H^ff^, the expressions of (ImA2)' and ImAo in terms of Wilson coeffi- 
cients and of the S-parameters are obtained 

Wo = -^lm{v*Vtd)[-(CeBe + ^C,B,^Z + {C4B^+ 

l^B,) X + CBy^ (f + f + f ) + CsBs^' (2>^+ 
l,^yC.,y^^,(9f,c.B,)x}, (71) 



(W2)' = -G^Im(v-,:V-,,) {CySf (I - I) + 

CsBr{Y-^)+C,Br'f]. (72) 

Notice that the matrix elements of the electromagnetic left-right operators 
Qt.S) which belong to the (8l, Sr) representation of SU (3) l'S>SU (3)^, contain 
Y and do not vanish in the chiral limit. 

The evaluation of the 5-factors requires a non-perturbative technique. 
The Wilson coefficients and the hadronic matrix elements both depend on 
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the regularization scheme. In order to cancel this dependence (up to 0(a^)), 
it is necessary to control the matching between the 5-parameters and the 
coefficients at the next-to-leading order. Notice that many non-perturbative 
methods (e.g. expansion) do not fulfil this requirement. 

Two different approaches to the matrix element evaluation have been 
used in recent next-to-leading e'/ e analyses: 



In our previous analysis [§, |] , the numerical values of the _B-parameters 



have been taken from lattice calculations Suitable renormalization 
factors are introduced to take into account the difference between the 
HV, NDR and lattice regularization schemes. For those 5-parameters 
not yet computed on the lattice^, we have made educated guesses, which 
are discussed in detail in ref. [Bl. 



A phenomenological approach has been implemented in ref. [0] , where 
the i?-parameters are constrained by using the experimental informa- 
tion from CP-conserving processes, by assuming SU{3) flavour symme- 
try and deducing some constraints relating hadronic matrix elements 
at the charm threshold. Unfortunately, there is no way to determine 
the most important 5-factors necessary to estimate e'/e, namely Bq 
and Bg, which remain essentially unconstrained in this approach. 



8 Results 



In this section, the main results of our analysis are summarized. These 
results have been obtained by varying the experimental quantities, e.g. the 
value of the top mass mt, tb, etc., and the theoretical parameters, e.g. the B- 
parameters, the strange quark mass m^, etc., according to their errors. Values 
and errors of the input quantities used in the following are reported in tables 
We assume a Gaussian distribution for the experimental quantities and 
a flat distribution (with a width of 2a) for the theoretical ones. The only 
exception is m^, taken from quenched lattice QCD calculations, for which we 



have assumed a Gaussian distribution, according to the results of ref. ||32 



^Indeed i?-parameters, which give the main contribution to the value of e'/e, namely 
Bq, B^^^^ and Bg^^^\ have already been computed on the lattice, see table ^. 
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Parameters 


Values 




(174 ± 17) GeV 


m,(2 GeV) 


(128 ± 18) MeV 


^^QCD 


(230 ± 80) MeV 


Kfc = AX' 


0.040 ± 0.006 


\Vub/Vcb\ = Ao- 


0.080 ±0.015 


tb 


(1.49 ±0.12) X 10-12 




0.685 ±0.076 


{fBBy^)th 


(200 ± 40) MeV 


VLib 


0.25 ±0.10 



Table 2: Values of the fluctuating parameters used in the numerical analysis. 



The theoretical predictions (cos 5, e'/e, etc.) depend on several fluctu- 
ating parameters. We have obtained their distributions numerically, from 
which we have calculated the central values and the errors reported below. 

Using the values given in the tables and the formulae given in the previous 
sections, we have obtained the following results: 

a) The distribution for cos 5, obtained by comparing the experimental 
value of e with its theoretical prediction, is given in fig. |^. As already 
noticed in refs. |]^, § and large values of and rrit favour 
cos 5 > 0, given the current measurement of Xd- When the condition 
160 MeV < fBB]l'^ < 240 MeV is imposed (/s-cut), most of the neg- 
ative solutions disappear, giving the dashed histogram of fig. 0, from 
which we estimate 

cos 5 = 0.47 ±0.32 . (73) 

b) A contour plot in the p-r] plane is given in fig. |[ It shows the current 
limits on the unitarity triangle defined in fig. |I|. 

c) In fig. ^, several pieces of information on e'/ e are provided. Lego plots 
of the distribution of the generated events in the e'/e-cos5 plane are 
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m,= (l74±17) GeV A'**ocd = (330 ± 1 00) MeV 
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Figure 3: Distributions of values for cos6, sin 2/3, e'/e and e' / e{A^ri)'^ , for 
rrit = (174 ± 17) GeV, using the values of the parameters given in tabs. 
The solid histograms are obtained without using the information coming from 
Bd-Bd mixing. The dashed ones use the Xd information, assuming that 160 
MeV< feB]^^ < 240 MeV. 
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m,= (l74±17) GeV A'**ocd = (330 ± 1 00) MeV 
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Figure 4: Contour plots in the plane. The solid, dashed and dotted 
contours contain 5%, 68% and 95% of the generated events respectively. The 
contours are given by excluding or including the fB-cut. Similar results can 
be found in refs. [2di , [^/ . 
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m,= (174±17) GeV A'*'ocd=(330± 1 00) MeV 
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Figure 5: Distributions of the events in the plane e'/e-cos5 without and with 
the fB-cut. The corresponding contour plots are displayed below the Lego 
plots. 
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Constants 


Values 


Gf 


1.16634 X IQ-^GeV"^ 


ITLc 


1.5 GeV 




4.5 GeV 


Mw 


80.6 GeV 




140 MeV 


Mk 


490 MeV 


Mb 


5.278 GeV 


AMk 


3.521 X 10-^2 MeV 


U 


132 MeV 


h 


160 MeV 


X = sin 9c 


0.221 


^exp 


2.268 X 10-3 


ReAo 


2.7 X 10-^ GeV 




0.045 




2 GeV 



Table 3: Constants used in the numerical analysis. 



shown, without and with the fs-cut. In the same figure, the corre- 
sponding contour plots are displayed. One notices a very mild depen- 
dence of e'/e on cos 5. As a consequence, one obtains approximately 
the same prediction in the two cases (see also fig. |^). In the HV scheme 
the results are 

e'/e = (2.3 ± 2.1) x 10"^ no - cut (74) 

and 

e'/e = (2.8 ± 2.4) x 10"^ fs - cut , (75) 
whereas in the NDR scheme we obtain 

e'/e = (2.8 ± 2.2) x 10"^ no - cut (76) 

and 

e'/e = (3.4 ± 2.5) x 10"^ /b - cut . (77) 
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^^3/2) 


Bl-2 


-83,4 


-85,6 


r(1/2) 


^7-8 


0.75 ±0.15 


0.62 ±0.10 


- 0.15^*) 


1 - 6^*^ 


1.0 ±0.2 


iW 


1.0 ±0.2 



Table 4: Values of the B -parameters, for operators renormalized at the scale 
fi = 2 GeV. The only exception is B^j^\ which is the renormalization group 
invariant B -parameter. Bg has been taken equal to Bk, at any renormaliza- 
tion scale. The value reported in the table is B^^'^^jj, = 2 GeV). Entries with 
a are educated guesses, the others are taken from lattice QCD calculations. 



By averaging the results given in eqs. (ffq ) and (O); obtain our best 
estimate 

e'/e = (3.1 ± 2.5 ± 0.3) x 10"^ - cut , (78) 

where the third error comes from the difference of the central values 
in the two schemes and gives an estimate of the uncertainty due to 
higher-order corrections. 



A similar resuh has been obtained in ref. [Q, using a different approach 
to the hadronic-matrix-element evaluation. They quote 



e 



e 



(6.7 ± 2.6) X 10"^ (79) 



for mt = 130 GeV. For this value of the top mass, the cancellation between 
penguin and electropenguin contributions is less effective, thus their e'/ e pre- 
diction is significantly larger than ours. Actually the two predictions agree, 
once the difference in the top mass is taken into account^. It is reassuring 
that theoretical predictions, obtained by using quite different approaches to 
matrix elements evaluation, are in good agreement. 

On the basis of the latest analyses, it seems very difficult that e'/e is 
larger than 10 x 10~^. Theoretically, this may happen by taking the matrix 
elements of the dominant operators, and Qs, much more different than it 
is usually assumed. One possibility, discussed in ref. |T1], is to take i?6 ~ 2 



and -Bg ~ 1, instead of the usual values .Bg ~ -Bg ~ 1- To our knowledge, no 
coherent theoretical approach can accommodate such large values of Bq. 

^Our analysis gives e'/e = (6.3 ± 2.3) x 10"-* for mt = 130 GeV. 
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